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Abstract 

Background: Efficient production of SA in Escherichia coli has been achieved by modifying key genes of the central 
carbon metabolism and SA pathway, resulting in overproducing strains grown in batch- or fed-batch-fermentor 
cultures using a complex broth including glucose and YE. In this study, we performed a GTA to identify those genes 
significantly upregulated in an engineered £ coli strain, PB12.SA22, in mid EXP (5 h), early STA (STA1, 9 h), and late 
STA (STA2, 44 h) phases, grown in complex fermentation broth in batch-fermentor cultures. 

Results: Growth of £ coli PB12.SA22 in complex fermentation broth for SA production resulted in an EXP growth 
during the first 9 h of cultivation depending of supernatant available aromatic amino acids provided by YE because, 
when tryptophan was totally consumed, cells entered into a second, low-growth phase (even in the presence of 
glucose) until 26 h of cultivation. At this point, glucose was completely consumed but SA production continued 
until the end of the fermentation (50 h) achieving the highest accumulation (7.63 g/L of SA). GTA between EXP/ 
STA1, EXP/STA2 and STA1/STA2 comparisons showed no significant differences in the regulation of genes encoding 
enzymes of central carbon metabolism as in SA pathway, but those genes encoding enzymes involved in sugar, 
amino acid, nucleotide/nucleoside, iron and sulfur transport; amino acid catabolism and biosynthesis; nucleotide/ 
nucleoside salvage; acid stress response and modification of IM and OM were upregulated between comparisons. 

Conclusions: GTA during SA production in batch-fermentor cultures of strain PB12.SA22 grown in complex 
fermentation broth during the EXP, STA1 and STA2 phases was studied. Significantly, upregulated genes during 
the EXP and STA1 phases were associated with transport, amino acid catabolism, biosynthesis, and nucleotide/ 
nucleoside salvage. In STA2, upregulation of genes encoding transporters and enzymes involved in the synthesis 
and catabolism of Arg suggests that this amino acid could have a key role in the fuelling of carbon toward SA 
synthesis, whereas upregulation of genes involved in pH stress response, such as membrane modifications, suggests 
a possible response to environmental conditions imposed on the cell at the end of the fermentation. 

Keywords: Shikimic acid production, Escherichia coli PTS- strain, Batch-fermentor culture, Complex fermentation 
broth, Global transcriptomic analysis, Microarrays, Regulatory network 




MICROBIAL CELL 
FACTORIES 



* Correspondence: adelfo@ibt.unam.mx 

^epartamento de Ingeniena Celular y Biocatalisis, Institute de Biotecnologia, 
Universidad Nacional Autonoma de Mexico, Av. Universidad 2001, Col. 
Chamilpa, Cuernavaca, Morelos 62210, Mexico 
Full list of author information is available at the end of the article 

O© 2014 Cortes-Tolalpa et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the 
BiolVlGCl C6ntT8l Creative Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, 
distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public 
Domain Dedication waiver (http://creativecommons.Org/publicdomain/zero/1.0/) applies to the data made available in this 
article, unless otherwise stated. 



Cortes-Tolalpa et al. Microbiol Cell Factories 2014, 13:28 
http://www.microbialcellfactories.eom/content/13/1/28 



Page 2 of 20 



Background 

The SA pathway is the common route leading to the 
biosynthesis of aromatic compounds in bacteria and in 
several eukaryotic organisms such as ascomycetes fungi, 
Apicomplexa, and plants [1,2]. In Escherichia coli, the 
first step in this pathway is the condensation of the 
CCM intermediates PEP and E4P, into DAHP by the 
DAHP synthase isoenzymes AroF, AroG, and AroH, 
which are encoded by the aroF, aroG, and aroH genes, 
respectively (Figure 1). The DHQ synthase, encoded by 
aroB, converts DAHP into DHQ. Subsequently, DHQ 
dehydratase, encoded by aroD, converts DHQ into DHS, 
and this compound is then transformed to SA by shi- 
kimate dehydrogenase, which is encoded by aroE. Shi- 
kimate kinases I and II, encoded by aroL and aroK, 
respectively, convert SA into S3P. Finally, S3P is con- 
verted to EPSP by 3-phosphoshikimate-l-carboxyvinyl- 
transferase, which is encoded by aroA. The last step in 
the SA pathway is the synthesis of CHA by the CHA 
synthase enzyme, which is encoded by aroC. CHA is the 
common building block for the formation of aromatic 
amino acids and compounds such as quinone, mena- 
quinone and enterobactin [3-5]. 

SA is a commercially important compound because it 
is considered to be an enantiomerically pure building 
block that is used as the precursor for the synthesis of 
numerous chemicals. Currently, SA has gained great im- 
portance as the starting compound for the chemical syn- 
thesis of OSP, the selective and potent inhibitor of the 
neuraminidase enzyme located on the surface of the 
influenza virus, known commercially as Tamiflu® and 



produced by Roche Pharmaceuticals [3-5]. OSP prevents 
the release of newly formed virus particles from influ- 
enza virus types A and B, avian influenza virus H5N1 
and, recently, human influenza virus H1N1. Since 1999, 
Roche Pharmaceuticals increased the production of OSP 
to ensure a significant reservoir in several countries in 
anticipation of a possible pandemic influenza outbreak; 
however, it has been estimated that in this scenario, the 
production of the antiviral would be insufficient to cover 
the requirements of the world population [6,7], particu- 
larly in developing countries such as Mexico. The latest 
human influenza outbreak, which appeared in Mexico in 
2009, showed that production of OSP is clearly insuffi- 
cient to satisfy the demand for this antiviral in an emer- 
gency situation. Additionally, the main supply of SA for 
OSP production is currently derived from the seeds of 
Chinese star anise (Illicium verum). The supply of this 
source is susceptible to vagaries of the weather. The star 
anise plant takes around six-years from planting to bear 
fruit but remains productive for a long time; additionally, 
extraction and purification from its seeds results in low 
yields. Thus, alternative biotechnological strategies with 
engineered bacterial strains to produce SA have gained 
relevance [4,5]. 

Metabolically engineered E. coli strains for SA produc- 
tion include several genetic modifications in CCM such 
as the introduction of an additional plasmid-copy 
DAHPS AroF^, encoded by arof° T or DAHPS AroG^ 
(aroG^), the tktA gene encoding transketolase I, and 
genes encoding enzymes from the SA pathway, including 
the single or double inactivation of genes aroK and aroh, 
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Figure 1 Central carbon metabolism and shikimic acid pathways in evolved E. coli PB12.SA22 strain. Glucose transport and 
phosphorylation are performed by galactose permease (GalP) and glucokinase (Glk), respectively. Genes and coded enzymes: golP, GalP; glk, Glk; 
tktA, transketolase I; pykF, pyruvate kinase I; pykA, pyruvate kinase II; ppsA, phosphoenolpyruvate synthase; aroF, aroG, aroH, DAHP synthase 
isoenzymes F, G and H, respectively; aroB, DHQ synthase; aroD, DHQ dehydratase; aroE, shikimate dehydrogenase; aroK, shikimate kinase I; aroL, 
shikimate kinase II; oroA, EPSP synthase; oroC, chorismate synthase. p}LBoroG fbr tktA and pJOPOoroBoroD, indicate cloned plasmids for SA 
overproduction in this strain [10]. Continuous arrows represent unique reactions catalyzed by one or more enzymes; dotted lines or arrows 
represent two or more enzymatic reactions or incomplete characterized reactions. Cross line in the reaction of SHK to S3P indicates interruption 
of the SA pathway. 
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and the introduction of an additional plasmid-copy of 
genes encoding limiting enzymes of the pathway such as 
aroB and aroE, resulting in an increased carbon flux 
from the CCM intermediates PEP and E4P to the SA 
pathway and accumulation of SA. The above-described 
genetic modifications in specific E. coli genetic back- 
grounds with additional modifications and grown under 
diverse culture conditions have resulted in the successful 
overproduction of SA with yields ranging from 0.08 to 
0.42 mol SA/mol glucose [8-13] (Table 1). 

Previous characterization of SA production in strain 
PB12.SA22 in 0.5-L batch-fermentor cultures using com- 
plex fermentation broth including 25 g/L of glucose and 
15 g/L of YE showed a characteristic two-phase growth 
behavior with an initially EXP growth with high [i while 
consuming ~ one-third of the initially added glucose and 
low level production of SA and other pathway interme- 
diates. During the second growth phase, the \a decreased, 
and the culture entered the STA phase despite the pres- 
ence of abundant residual glucose in the supernatant 
broth, whereas SHK pathway intermediate production 
increased, continuously reaching its maximum until the 
end of the fermentation (50 h). Interestingly, residual 
glucose was depleted from supernatant culture during 
the STA phase associated with SA production [10]. This 
growth, glucose consumption and SA production behav- 
ior suggest that during the EXP growth phase, strain 
PB12.SA uses YE components to support growth and, as 
consequence of the possible depletion of an essential nu- 
trient component, cell ceases growth upon entering the 
STA phase, where residual glucose was channeled by this 
strain to produce SA and other pathway intermediates. 
These data suggest the presence of important genetic 
regulation and physiological differences during the EXP 
and STA phases. 



GTA has been proven to be a powerful tool to study 
regulation of cellular metabolism in response to specific 
environmental conditions. To our knowledge, in relation 
to SA production in E. coli, GTA has been used to eluci- 
date byproduct formation in the SA production strain 
W3110.shikl under carbon and phosphate-limited (car- 
bon-rich) chemostat conditions, suggesting that bypro- 
duct formation under carbon limitation is explained by 
the upregulation of a set of genes coupled to the SA 
pathway [14]. In this study, we report the GTA in strain 
PB12.SA22 during SA production in batch-fermentor 
cultures using complex fermentation broth including 
25 g/L of glucose and 15 of g/L YE. GTA was deter- 
mined in this strain in the mid EXP growth phase (5 h), 
at the beginning of the STA phase (STA1) (9 h) and in 
the late STA phase (STA2) (44 h), and the comparison 
of differentially upregulated genes was established during 
the EXP/STA1, EXP/STA2 and STA1/STA2 phases to 
correlate changes in the global expression profile and 
growth, glucose consumption and SA production pro- 
files in this strain. 

Results 

Growth and SA pathway intermediate production 

Based on growth and glucose consumption profiles, 
strain PB12.SA22 showed its characteristic two-phase 
growth behavior [10], with an initial EXP growth phase 
during the first 8 h of cultivation reaching an OD 600 nm = 
13.37 with a t* = 0.47 ±0.002 h" 1 and a q s = 334± 
0.29 mmol glucose g DW" 1 h" 1 . After 8 h of cultivation, 
the culture reduced its growth rate even in the presence 
of -70% of initially added glucose in the fermentation 
broth, and after 10 h of cultivation, the strain entered 
into a pseudo-STA phase in which a low growth phase 
was observed until 26 h of cultivation, where the highest 



Table 1 E. coli SA engineered overproducing strains, growth conditions and SA yield 


Strain/derivative 


Relevant characteristics 


Culture conditions 


SA 

[g/L] 


Yield (mol 
SA/mol glucose) 


Reference 


SP1.1/pKD12.138 


serA-aroB AaroL AaroK pSV]8aro? br P tac aroE 
serA tktA 


1-L fed-batch cultures, mineral broth 
with 25 g/L of glucose and 15 g/L of YE 


52 


0.18 


[8] 


SP1.1pfs/pSC6.090B 


PTS" serA::oroB AoroL AaroK PTS" P tac glf glk 
aro? br tktA P tac aroE serA, glf 1 glk 2 


10-L fed-batch reactors, mineral broth 
with 25 g/L of glucose and 15 g/L of YE 


71 


0.27 


[8] 


W31 lO.shikl 


AaroL aroG fbr aro^tnaA, and plasmid 
overexpressed aroK 


Chemostat cultures using mineral broth 
under glucose limiting conditions 


NR 


-0.08 


[9] 


PB12.SA22 (JM101 
derivative) 


PTS" AaroL AaroK pJLBaroG fbr tktA pTOPO 
aroB aroE 


0.5-L batch reactors, mineral broth with 
25 g/L of glucose and 1 5 g/L of YE 


7 


0.29 


[10] 


DHPYAAS-T7(DH5a- 
derivative) 


Aptstilcrr AaroL Aarok AydiB knock-in of 
T7-RNA-pol gene, pAOC-TGEFB aroE aroB, 
site-specific mutagenesis glk tktA arof° r 


Fed-batch fermentation, modified M9 
medium, with 25 g/L of glycerol and 
25 g/L of YE 


1.85 


NR 


[11] 


AR36 (JM101 
derivative) 


PB12 lacl~ a rotC aroL~ pykF Trc/aroB + 
tktA + aroG fbr+ aroE + aroD + zwf + 


Batch fermentation, mineral broth with 
1 00 g/L of glucose and 30 g/L of YE 


43 


0.42 


[12] 


E. coli SA1 16 


Chromosomally evolved and cofactor 
metabolic engineered strain 


Mineral broth with 10 g/L of glucose, 
1 g/L of peptone, 1 .24 g/L of proline 


3.12 


0.33 


[13] 



Glucose facilitator {glf) and 2 Glucokinase (glk) from Zymomonas mobilis; NR, Non-reported. 
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OD 60 o nm = 15.77 was observed. During this period, the 
strain completely consumed the residual glucose present 
in the supernatant broth (showing a q s = 2.68 ± 0.5 mmol 
glucose g DW" 1 h" 1 ). From this point to the end of the 
fermentation, the biomass showed a slight reduction to 
OD 600 nm ~ 15 (Figure 2A). Interestingly, analysis of the 
consumption profile of available supernatant aromatic 
amino acids by HPLC showed that strain PB12.SA22 to- 
tally consumed the tryptophan present in the YE during 
EXP growth, reducing its growth when this aromatic 
amino acid was completely depleted from the culture 
supernatant (Figure 2B). SA production was detected in 




0 10 20 30 40 50 60 
Time (h) 




Time (h) 

Figure 2 Growth, aromatic intermediates production, and 
substrate consumption during batch fermentation cultures of 
E. coli PB12.SA22 strain grown in complex broth. (A) Biomass 

production (•) Glucose consumption (♦). (B) Residual Phenylalanine 
(□), Tyrosine (o), Tryptophan (a). (Q SA (A), DHS (♦), GA (T) 

production. In panel A, numbers in parenthesis indicate the 
sampling time of biomass for total RNA extraction used for 
microarray analysis (1), (2) and (3) indicates samples collected at 5, 9, 
and 44 h of cultivation, respectively. 



supernatant cultures in the EXP phase, but the highest 
SA accumulation was observed during the last 24 h, ac- 
cumulating 7.63 g/L SA at the end of fermentation 
(50 h) (Figure 2C). The determined SA pathway bypro- 
ducts, DAHP, DHS, and GA, appeared to be associated 
with growth because their production remained rela- 
tively constant after strain PB12.SA22 entered the STA 
phase. SA production resulted in a yield of 0.31 mol SA/ 
mol glucose, and the total aromatic yield (the combined 
yield of SA, DAHP, DHS, and GA) was 0.36 mol aro- 
matic intermediates/mol glucose. 

Global transcriptomic analysis during the EXP and STA 
phases 

GTA analysis was determined by triplicate in the middle 
of the EXP phase (5 h), at the end of the EXP growth 
phase and upon entering pseudo-STA phase (9 h) 
(STA1) and at the end of the fermentation, or the late 
STA phase (44 h) (STA2). Average expression data were 
compared among EXP/STA1, EXP/STA2, and STA1/ 
STA2, using the expression data of 4070 genes from the 
array of E. coli MG1655 included in the Affymetrix Gen- 
eChip® E. coli Genome 2.0 by the RP method. For all ex- 
pression data, the RP method calculated a FDR. Those 
genes with an FDR value = 0 have the highest probability 
of biological relevance [15,16]. To identify differentially 
expressed (upregulated and downregulated) genes among 
compared growth conditions, an FDR value cutoff < 0.05 
was used. This analysis resulted in the identification of the 
significant upregulation of 43 genes and the downregula- 
tion of 80 genes in EXP/STA1, the upregulation of 37 
genes and the downregulation of 38 genes in EXP/STA2, 
and the upregulation of 50 genes and the downregulation 
of 47 genes in ST1/STA2 (Table 2). Among them, 22 com- 
mon genes were found to be upregulated during the entire 
fermentation process (EXP, STA1 and STA2). Biological 
functions were assigned to upregulated genes according to 
the EcoCyc database for E. coli strain MG1655 and 
broadly grouped on the basis of GO terms [17]. 

Initial GTA data mining for upregulated genes in 
EXP/STA1 (Table 3), EXP/STA2 (Table 4), and STA1/ 
STA2 (Table 5) comparisons resulted from the RP ana- 
lysis from average transcription data for genes encoding 
enzymes involved in glycolysis, PPP, TCA, and the glyox- 
ylate shunt, acetate metabolism and gluconeogenic cap- 
abilities in the EXP, STA1 and STA2 phases, showed that 
genes encoding proteins of the CCM and SA pathways 
were not significantly upregulated among performed com- 
parisons (Additional file 1). However, significantly upregu- 
lated genes detected during the EXP and STA phases 
suggest important differences in sugar, amino acid, nu- 
cleotide/nucleoside and ion transport; metabolic process, 
particularly amino acid catabolism and biosynthesis; nu- 
cleotide/nucleoside salvage; acid stress response and 
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Table 2 Differentially expressed genes during EXP/STA1, EXP/STA2 and STA1/STA2 comparisons during SA production 9 

GOterms b EXP/STA1 EXP/STA2 STA1/STA2 





UP 


DR 


UP 


DR 


UP 


DR 


Transport 


15 


9 


9 


5 


8 


10 


Amino acid metabolism 


9 


15 


8 


ND 


16 


4 


Electron carrier activity 


2 


ND 


1 


ND 


3 


ND 


Out of GOs 


2 


21 


4 


8 


4 


14 


ATP catabolic/biosynthetic process 


2 


1 


ND 


ND 


ND 


ND 


Nucleotide/nucleoside metabolism 


6 


2 


1 


ND 


ND 


2 


Sulfur metabolism 


2 


ND 


4 


ND 


2 


ND 


Stress response 


1 


12 


2 


6 


5 


3 


Catalytic activity 


1 


ND 


1 


ND 


ND 


ND 


Regulatory 


2 


7 


4 


7 


4 


4 


Cell envelope 


ND 


5 


2 


7 


1 


4 


Structural molecule activity 


1 


ND 


1 


ND 


ND 


ND 


Lipid metabolism 


ND 


1 


ND 


ND 


3 


ND 


Other metabolic process 


ND 


7 


ND 


5 


4 


6 


Total gene number 


43 


80 


37 


38 


50 


47 



a Based on resultant FDR value from the RP analysis equal or less to 0.05 [15]; b Biological functions were assigned according to the EcoCyc database for £ coli 
strain MG1655 and broadly grouped on the basis of GO terms [17]. ND, Non-detected. 



modification of the IM. A detailed list of those differen- 
tially expressed genes in the EXP/STA1, EXP/STA2 and 
STA1/STA2 comparisons, expression data and FDR 
values information are presented in Additional file 1. 
The role and the relationship with growth, substrate 
consumption and SA production of relevant upregu- 
lated genes during EXP/STA1, EXP/STA2 and STA1/ 
STA2 will be discussed below. 

Discussion 

Upregulation of genes involved in sugar transport 

Genes encoding sugar transporters that were found to 
be upregulated in the EXP/STA1 comparison included 
the sugar porin LamB (lamB), the periplasmic binding 
protein, the membrane subunit and the ATP binding 
subunit components of the ABC maltose transport sys- 
tem (malE and rnalFK, respectively) [17] and the GlpT 
glycerol 3-P MFS transporter (glpT), whereas in the 
EXP/STA2 comparison, the upregulation of the non- 
specific OMP ompC (OmpC), a general porin and the 
periplasmic binding protein of the maltose ABC 
transporter (malE) was observed. A previous transcrip- 
tomic analysis of our group determined by RT-qPCR in 
cultures of the parental strain PB12 in M9 minimal 
broth showed the upregulation of the galP gene, suggest- 
ing that GalP protein is the main transporter used for 
glucose import in this PTS" strain [18], as is the case in 
a derived PB12 strain for L-Phe production grown in 
M9 broth supplemented with 5 g/L of YE [19]. However, 



our GTA showed that in the SA-producing strain PB12. 
SA22, the galP gene was not significantly upregulated 
during EXP growth and the STA phase. This result was 
also observed in GTA for the parental strain PB12 in the 
same culture conditions (data not shown), suggesting 
that in both strains grown in complex fermentation 
broth for SA production, the OMP OmpC, LamB, and 
the maltose ABC transporter, have an important role in 
the transport of glucose because previous works re- 
ported the participation of these proteins in the uptake 
of glucose in E. coli [20-23]. All these genes were found 
to be upregulated in the EXP and STA1 phases, associated 
with the consumption of >30% of the initially added glu- 
cose to the fermentation broth; however, only in the ST All 
STA2 comparison was detected the upregulation of ompC. 
Differential expression of these transporters is in agree- 
ment with the absence of glucose in the supernatant 
culture observed at the middle STA phase (26 h). Interest- 
ingly, the gene encoding the glycerol-3-P (G3P) MFS trans- 
porter (GlpT) (glpT) was also detected upregulated during 
the EXP/STA1 comparison but not during STA1/STA2. 
G3P plays a major role in glycolysis and phospholipid bio- 
synthesis in E. coli. G3P is transported by GlpT and is 
reduced by aerobic or anaerobic G3P dehydrogenase 
into dihydroxyacetone phosphate, which is converted into 
fructose- 1,6-diphosphate or glyceraldehyde-3-phosphate 
and then enters the glycolysis pathway. In phospholipid 
biosynthesis, G3P forms the backbone of all phospholipid 
molecules and the polar groups of phosphatidylglycerol 
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Table 3 Average expression data of upregulated genes during EXP/STA1 comparison based on resultant FDR value 
from the RP analysis equal or less to 0.05 in compared growth phases during SA production cultures 



Gene 


Locus 


Protein name 3 


Cellular function 13 


Average 
expression 




b3744 


A^nprpninp Qvnthpfp^p A 
AAbpdldyil Ic by 1 III Icldbc AA 


Aminn priH mpfphnlKm 

/\l 1 III IU CIL.IU 1 1 ICluUUIIjl 1 1 


28 5470 


asp A 


hzll 


Abpdlldlc dllllllUMId lydbc 


AllllllU dLIU II IcldUUIIbl 1 1 


J . i *+*+o 




b2904 


(nk/rinp rlp^x/^np w^fpm H nrnfpin 
vjiy^_n ic Licavayc jyoici i i i i uiuicii i 


Aminn ^riH mpt^hnli^m 
i\\ \ in \\j aLiu 1 1 ic la u\j\ ioi 1 1 


7.2457 


Q gcvP 


b2903 


f^k/rinp Hprprhnwlp^p 

OiyLIMc UcLal UUAy Idbc 


Aminn priH mpfphnlKm 
/\i i in \kj ciL-iu 1 1 iciauuiiji 1 1 


1 2 2509 


gcvT 


b2905 


Aminnmpthvltr^n^fpr^^p 

l\\ 1 III IUI 1 IC 11 lyiLIQI 1 0 1 CI COC 


Amino ?iriH mpt3hnli<;m 
r\\ \ hi iu 0\^\\J i i iciakjwiioi i i 


1 2.0966 


put A 


b1014 


Fused PutA DNA-binding transcriptional repressor/proline dehydrogenase/ 

I py 1 I Ul II Ic J Ldl UUXy Id Lc Ucl lyUIUycl Idbc 


Amino acid metabolism 


11.8149 


sdoB 


b2797 


I -cprinp Hppminpc:p II 

l_ jCI II IC UCCII I III IOjC II 


Aminn priH mptphnli^m 

/\l 1 III IU CIL.IU 1 1 ICICIUUIIjI 1 1 


8 2308 


tnnA 
Li IUr\ 


Uj / uo 


L Lyblclllc UcbUII liyUldbc/ LiypLUpi ldl Idbc 


AllllllU dLIU 1 1 IcldUUIIbl 1 1 




LI lUL 


h3707 
Uj 1 \J 1 


LI IU UptrlUII IcdUtrl ptrpilUc 


AllllllU dLIU 1 1 IcldUUIIbl 1 1 


91 7140 


QtpG 


Uj/jj 


ATP c\ /r\1~h\ 2 ca P r~r\ rv~\ /^i^rv^rv^^ ci il^i irul - 
Air by 1 1 LI Idbc r-| LUl 1 ipicX ydlllllld bUUUl 11 1 


Air UlUbyi III IcllL/LdldUUIIL piULcbb 




fecE 


b4287 


rclllL UlLllldlc nDL UdllbpUllcl Air UlllUlliy bUUUIIIl 


ATP Y\\c\z\ir\XY\c±X\c 1 o\~> Kol \c pitpippcc 
Air UlUby 1 1 11 IcllL/ Ld IdUUI IL piULcbb 


7 01 84 


ompT 


UUOOJ 


OK/1 nrntaaca \/ll (f^\!\A nrr^toin 
wlvl piULcdSc VII ^wlvl piULclll Du) 


C \\ iX\ r d r~X\\ /'\X\ i 
L.dldiyilL dLllvliy 


O.UOp'O 


groS 


B4142 


rUiypcpLIUc. OIUHO, LlldpclUllc UIMUb LU nbpOU III picb. IViy Air, bUppicbblliy lib 

ATPase activity 


Cell processes 


5.9806 


cyoB 


b0431 


Cytochrome bo terminal oxidase subunit 1 


Electron carrier activity 


6.3920 


cyoD 


b0429 


Cytochrome bo terminal oxidase subunit IV 


Electron carrier activity 


6.6316 


gpt 


b0238 


Xanthine-guanine phosphoribosyltransferase 


Nucleotide/nucleoside metabolism 


6.4026 


guaC 


b0104 


GMP reductase 


Nucleotide/nucleoside metabolism 


6.7694 


pyrB 


b4245 


Aspartate carbamoyltransferase, catalytic subunit 


Nucleotide/nucleoside metabolism 


10.8487 


rihA 


b0651 


Ribonucleoside hydrolase 1 (pyrimidine-specific) 


Nucleotide/nucleoside metabolism 


5.9568 


upp 


b2498 


Uracil phosphoribosyltransferase 


Nucleotide/nucleoside metabolism 


5.4210 


yliF 


b0834 


Predicted diguanylate cyclase 


Nucleotide/nucleoside metabolism 


5.6595 


nlpA 


b3661 


Lipoprotein-28 


Out of GOs 


7.2359 


yiel 


b3716 


Predicted IM protein 


Out of GOs 


6.6363 


gcvB 


b4443 


GcvB small regulatory RNA 


Regulatory 


7.4984 


yeeN 


b1983 


Conserved protein 


Regulatory 


5.3142 


rihA 


b4142 


GroES, chaperone binds to Hsp60 in pres. Mg-ATP, suppressing its ATPase activity 


Stress response 


5.9806 


cyoE 


b0428 


Heme 0 synthase 


Structural molecule activity 


5.9183 


cysC 


B2750 


Adenylylsulfate kinase 


Sulfur metabolism 


7.0141 


cysH 


b2762 


3'-Phospho-adenylylsulfate reductase 


Sulfur metabolism 


11.2772 


cirA 


b2155 


0M receptor involved in uptake of ferric dihyroxybenzoylserine 


Transport 


9.4700 


fecA 


b4291 


Ferric citrate 0MP FecA 


Transport 


19.9369 


fecB 


b4290 


Ferric dicitrate ABC transporter - periplasmic binding protein 


Transport 


14.5864 


fruB 


b2169 


Fructose PTS permease - FruB subunit 


Transport 


5.8528 


glpT 


b2240 


GlpT glycerol-3-P MFS transporter 


Transport 


8.6190 


lamB 


b4036 


Phage lambda receptor protein; maltose high-affinity receptor 


Transport 


18.0411 


livK 


b3458 


Leucine ABC transporter - periplasmic binding protein 


Transport 


6.0726 


molE 


b4034 


Maltose ABC transporter - periplasmic binding protein 


Transport 


174.0512 


malF 


b4033 


Maltose ABC transporter - membrane subunit 


Transport 


19.1719 


molK 


b4035 


Maltose ABC transporter - ATP binding subunit 


Transport 


14.8026 


nmpC 


b0553 


0MP protein; locus of qsr prophage 


Transport 


11.7561 


proX 


b2679 


Glycine betaine/proline ABC transporter - periplasmic binding protein 


Transport 


7.3363 
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Table 3 Average expression data of upregulated genes during EXP/STA1 comparison based on resultant FDR value 
from the RP analysis equal or less to 0.05 in compared growth phases during SA production cultures (Continued) 

sbp b3917 Sulfate/thiosulfate ABC transporter - periplasmic binding protein Sbp Transport 1 1.2712 

tnaB b3709 TnaB tryptophan ArAAP transporter Transport 9.9672 

tsx b041 1 Nucleoside channel; receptor of phage T6 and colicin K Transport 6.2757 

Retrieved from EcoCyc database, b Biological functions were assigned according to the EcoCyc database for E. coli strain MG1655 and broadly grouped on the 
basis of GO terms [17]. 



and cardiolipin [24]. Upregulation of this gene in strain 
PB12.SA22 is intriguing because G3P was not present dur- 
ing the SA production cultures. GlpT could be involved 
in phospholipid biosynthesis during the EXP/STA1 
comparison because it was previously reported, and in- 
direct evidence suggests that GlpT expression may be 
regulated by a byproduct from the glycerolipid biosyn- 
thetic pathway [24]. 

Upregulation of genes involved in the amino acid 
transport and metabolism 

During the EXP/STA1 comparison, UvK, encoding a 
subunit of the ABC transporter for leucine, and tnaB 
(part of the tnaCAB operon) encoding the TnaB trypto- 
phan ArAAP transporter, which is proposed to be in- 
volved in tryptophan scavenging [25], were upregulated. 
The EXP/STA2 comparison revealed the upregulation of 
the genes for the putative membrane protein and the 
ATP binding components of the BppBCDF dipeptide 
ABC transporter (dppB and dppF), which is associated 
with the transport of proline-glycine as source of pro- 
line, histidine-glutamic acid, and leucine-tryptophan 
[26,27]; the putative membrane protein of the murein 
tripeptide ABC transport system (oppB), which can 
transport oligopeptides of up to five residues in length 
as well as recycle cell-wall peptides, but it has been re- 
ported to lack affinity for free amino acids [28-31]; and 
the complete glycine betaine/proline ABC transporter 
encoded by the proVWX operon [17]. Expression of this 
operon was reported to substantially increase at high 
osmolarity to scavenge glycine, betaine, proline, taurine, 
ectoine, carnitine as the precursor for betaine, and cho- 
line to achieve high intracellular concentrations of these 
osmoprotectants [32-34]. During the STA1/STA2 com- 
parison, the periplasmic binding protein of the arginine 
ABC transporter {art/); the 4-aminobutyrate antiporter 
glutamic acid GadC (gadC), a member of the APC 
superfamily of amino acid transporters [35]; and the in- 
tegral membrane subunit and the ATP-binding compo- 
nent of the lysine/arginine/ornithine ABC transporter 
and the histidine ABC transporter (hisM and hisP, re- 
spectively) were found to be upregulated [17]. 

Differential upregulation of those genes encoding pep- 
tide and amino acid transporters during the entire fer- 
mentation process suggests that during batch culture for 



SA production, strain PB12.SA22 imports the substrates 
that are available in the fermentation broth into the cell 
because they are supplied by YE, which contains 6% 
amino nitrogen, supplying peptides and 18 of the 20 
proteinogenic amino acids (with the exception of aspara- 
gine and glutamine) [36]; however, differential upregula- 
tion of amino acid transporters observed during the 
EXP, STA1 and STA2 phases suggests a possible differ- 
ential amino acid requirement by the cell during the 
fermentation process, particularly for leucine and trypto- 
phan during the EXP/STA1 comparison and leucine, ar- 
ginine, glutamic acid, lysine and histidine during the 
STA1/STA2 comparison. Amino acid availability in the 
fermentation broth provided by YE almost certainly alle- 
viates the cellular requirements because YE contains all 
these amino acids [36]. Finally, in the EXP/STA1 and 
EXP/STA2 comparisons, ornpT, encoding protease VII 
was upregulated; this protein is capable of cleaving sev- 
eral peptides at the center of paired basic residues but 
not at single basic residues, suggesting a distinct mech- 
anism for trypsin-like proteases [17]. Although OmpT is 
not involved in transport, it is possibly related to peptide 
assimilation in strain PB12.SA22 during SA production. 

Together with the upregulation of genes encoding 
peptide and amino acid transporters, several genes en- 
coding enzymes involved in amino acid catabolic or 
biosynthetic pathways were differentially upregulated 
during the EXP growth and STA phases. Among biosyn- 
thetic genes encoding the enzymes involved in the 
complete cysteine biosynthesis from sulfate, CysD, CysC, 
CysH (cysD, cysC and cysH, respectively) [17] were 
found upregulated during the EXP/STA2 and STA1/ 
STA2 comparisons, while CysJ {cisj), encoding sulfite 
reductase (NADPH), which catalyzes the 6-electron re- 
duction of sulfite to sulfide, one of several activities ne- 
cessary for the biosynthesis of cysteine from sulfate [17], 
was found to be upregulated only in the EXP/STA1 
comparison. During the EXP/STA2 comparison, we de- 
tected the upregulation of those genes encoding hydro- 
xymethyltransferase (glyA) involved in the conversion 
of serine to glycine; the enzyme cysteine synthase A 
(cysK) catalyzing the conversion of O-acetyl-L-serine to 
L-cysteine [37-39] (although, in the absence of a sulfur 
source, this enzyme catalyzes the slow conversion of O- 
acetyl-L-serine into pyruvate + acetate + ammonia, or the 
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Table 4 Average expression data of upregulated genes during EXP/STA2 comparison based on resultant FDR value 
from the RP analysis equal or less to 0.05 in compared growth phases during SA production cultures 



Gene 


Locus 


Protein name 3 


Cellular function G0 b 


Average 
expression 


ykfE 


b0220 


{ivy) Protein: inhibitor of vertebrate C-type lysozyme 


Amino acid metabolism 


29.2978 


cysK 


b2414 


Enzyme: cysteine synthase A 


Amino acid metabolism 


83.0796 


gcvH 


b2904 


Polypeptide: glycine cleavage system H protein 


Amino acid metabolism 


30.4955 


gcvP 


b2903 


Enzyme: glycine decarboxylase 


Amino acid metabolism 


33.3045 


gcvT 


b2905 


Enzyme: aminomethyltransferase 


Amino acid metabolism 


38.3509 


gdhA 


b1761 


Enzyme: glutamate dehydrogenase 


Amino acid metabolism 


20.9785 


glyA 


b2551 


Enzyme: serine hydroxymethyltransferase 


Amino acid metabolism 


23.7380 


tnoA 


_Uj/Uo 


L LySLcll Ic UcSUII I lyUldSc/ LI ypLUpi Idl IdSc 


AlIllllU dLIU 1 1 IcLdUUIISI 1 1 


1 J7.D 1 ZJ 


tnoL 


b3707 


LI IU U|JclUII IcdUcI |Jc|JLIUc 


AlIllllU dLIU 1 1 IcLdUUIIol 1 1 


UZ.7/ ou 




b0565 


Fn7\/mp - OM nrnfp^p \/ll (OM nrnfpin 

l_IIZ.yillC. WIVI pikJlCQoC VII [-JlkJLCIII UUJ 


f^t^k/tir 3rti\/it\/ 
v-a laiy u^. a^-iiviiy 


35.7201 


fimC 


b4316 


PnlvnpntiHp - nprinl^mir rh?\nprrir\p rpm lirprl for tvnp 1 fimhri^p 
r \j\y pep Liuc. pci i \j\a o\ i 1 1^. ^.i lupci , u\ ic, i crunch Lypc i iniikjiiac 


f~p|| pn\/p|nnp 
^.cn ci ivci^jpc 


63.3469 


fimD 


b431 7 


PnK/npntiHp" OM nrnfpin - pynnrt ^nH accpmhlu nf tvnp 1 fimhri^p 
ruiypcpLiuc. wivi piuicni, ca|jui l ai iu ujjci i ikjiy ui Lypc i iniikjiiac 


f~p|| pn\/plnnp 
^.Cll ci ivciwpc 


25.6637 


nlpA 


b3661 


Pnk/npntiHp" linnnrnfpin-7R 
ruiypcpiiuc. 1 1 pupi w lci i i z.o 


f~p|| pn\/plnnp 
^.cn ci ivciwpc 


48.31 94 


\/rl!D 


U I Oy/ 


rUlypcpLIUc. pULdLIVt: bUUUNIL Ul TUIvj TUIn I IdVUpi ULcll I 


CIcLLlUN Ldlllcl dLLIVILy 


Ay.yyjyJ 


nrdA 


b2234 


Protein: ribonucleoside diphosphate reductase 1, a subunit dimer 


IMULIcULlUc/ 1 IULIcUoIUc 

metabolism 


37.6063 


yciW 


b1287 


Polypeptide: predicted oxidoreductasePolypeptide: predicted oxidoreductase 


Out of GOs 


116.1438 


yiel 


b3716 


{cbrB) Polypeptide: predicted IM protein 


Out of GOs 


31.9090 


yeeN 


b1983 


Polypeptide: conserved protein 


Regulatory 


22.0828 


yhiE 


b3512 


{gadE} Polypeptide: GadE DNA-binding transcriptional activator 


Regulatory 


46.2077 


yhiW 


b3515 


(gadW) Polypeptide: GadW DNA-binding transcriptional dual regulator 


Regulatory 


37.7545 


hdeA 


b3510 


Protein: HdeA dimer, inactive form of acid-resistance protein 


Stress response 


32.4174 


hdeB 


b3509 


Polypeptide: acid stress chaperone 


Stress response 


56.2185 


b0023 


b0023 


Polypeptide: 30S ribosomal subunit protein S20 


Structural molecule activity 


45.8728 


cysC 


b2750 


Enzyme: adenylylsulfate kinase 


Sulfur metabolism 


52.8463 


cysD 


cysD 


Polypeptide: CysD 


Sulfur metabolism 


236.2573 


cysH 


b2762 


Enzyme: 3'-phospho-adenylylsulfate reductase 


Sulfur metabolism 


194.5029 


cysJ 




CllZyillc. SUIIILc IcUULLdSc, lldVUpiULcIN bUUUlllL LUIIipicX 


jUIIUI 1 IlcLdUUIIbl 1 1 


OZ.JJ^J 


ydjN 


b1 729 


PnlvnpntiHp - nrpHirfprl tr^n^nnrfpr 

r <U\y pcpLILJ C. piCLJI^-LCLJ LIUllop^lLCI 


Tr^n^nnrt 

1 1 Ul 1 J|JUI L 


43 4944 


dppB 


b3543 


Polypeptide: dipeptide ABC transporter - putative membrane subunit 


Transport 


23.0830 


dppF 


b3540 


Polypeptide: dipeptide ABC transporter - putative ABC binding subunit 


Transport 


32.6107 


fecA 


b4291 


Polypeptide: ferric citrate 0MP FecA 


Transport 


23.6130 


molE 


b4034 


Maltose ABC transporter - periplasmic binding protein 


Transport 


230.4651 


ompC 


b2215 


Transporter: 0MP C 


Transport 


61.2832 


oppB 


b1244 


Polypeptide: murein tripeptide ABC transporter/peptide ABC transporter - putative 
membrane subunit 


Transport 


44.3550 


proV 


b2677 


Polypeptide: glycine betaine/proline ABC transporter - ATP binding subunit 


Transport 


34.7488 


proW 


b2678 


Polypeptide: glycine betaine/proline ABC transporter - membrane subunitPolypeptide: 
glycine betaine/proline ABC transporter - membrane subunit 


Transport 


35.6683 


proX 


b2679 


Polypeptide: glycine betaine/proline ABC transporter - periplasmic binding protein 


Transport 


48.2957 



a Retrieved from EcoCyc database, b Biological functions were assigned according to the EcoCyc database for E. coli strain MG1655 and broadly grouped on the 
basis of GO terms [17]. 
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conversion of O-acetyl-L-serine into simple serine [40]; 
and glutamate dehydrogenase (gdhA) involved in the 
NADPH-dependent amination of a-ketoglutarate to 
yield L-glutamate. Catabolic genes upregulated during 
the EXP/STA1 and EXP/STA2 comparisons included 
tnaCA (part of the tnaCAB operon), which encodes the 
tryptophanase leader peptide and the tryptophanase in- 
volved in tryptophan degradation to indole + PYR + am- 
monia + H + . During the EXP/STA1 and EXP/STA2 
comparisons, the entire glycine cleavage system (GCV) 
(gcvTHP operon), a multienzyme complex catalyzing the 
reversible oxidation of glycine, yielding carbon dioxide, 
ammonia, 5,10-methylenetetrahydrofolate and a reduced 
pyridine nucleotide [17,41] was also upregulated. Tetra- 
hydrofolate serves as a recipient for Q units generated 
during glycine cleavage to form the methylene group, 
which are of central physiological importance, as will be 
discussed below [41-43]. 

Another group of upregulated genes found during the 
EXP/STA1 comparison, were the genes involved in am- 
monia dependent conversion of aspartate to asparagine 
(amino acid not supplied by YE), by asparagine synthe- 
tase (asnA) and the catabolic genes involved in the 
conversion of aspartate to Fum + ammonia + H + by 
aspartate ammonia-lyase {asp A), which is the first step 
in the proline degradation pathway to glutamic acid by 
proline dehydrogenase (putA), and in serine degradation 
to PYR + ammonia + H + by L-serine deaminase II (sdaB) 
[17]. Notably, a group of upregulated genes during the 
STA1/STA2 comparison include genes encoding the en- 
tire arginine biosynthesis pathway and the ornithine bio- 
synthesis pathway {argA, argCBH, argG, argD and argl) 
and carbamoyl phosphate synthetase (carAB), which cat- 
alyzes the first committed step in the biosynthetic path- 
ways for the production of arginine and pyrimidine 
nucleotides. The small subunit of this enzyme (carA) 
also hydrolyzes glutamine to glutamate and ammonia 
[44]. ilvG_l and UvGjL, part of the UvLXG_lG_2MEDA 
operon [45], are pseudogenes encoding subunits of acet- 
ohydroxybutanoate synthase / acetolactate synthase, an 
essential enzyme that catalyzes the biosynthesis of a- 
aceto-a-hydroxybutyrate for the isoleucine pathway and 
of a-acetolactate for the valine biosynthesis when E. coli 
K-12 was grown on acetate or oleate as the sole carbon 
source [46]. Other genes upregulated in the STA1/STA2 
comparison encode cobalamin-independent homocysteine 
transmethylase (metE), which catalyzes the final step of 
de novo methionine biosynthesis in the absence of ex- 
ogenously supplied vitamin B12 (cobalamin); a-keto 
glutarate reductase/D-3-phosphoglycerate dehydrogen- 
ase (serA), which catalyzes the first committed step in 
the biosynthesis of L-serine from 3-P-glycerate; and glu- 
taminase (ybaS), which catalyzes the degradation of glu- 
tamine to yield glutamate [17]. 



The main cellular role of amino acids is as structural 
blocks for the synthesis of proteins. Their availability de- 
termines the growth capabilities of cells, thereby defin- 
ing the extent of the EXP growth phase, particularly for 
auxotrophic strains. Because strain PB12.SA22 is auxo- 
trophic for aromatic amino acids due to the deletion of 
the aroK and aroL genes, its growth depends of the 
extracellular availability of these essential amino acids 
present in YE. As it was observed in Figure 2B, when 
tryptophan was depleted from the fermentation broth, 
this strain decreased its (A even in the presence of <70% 
residual glucose, showing that the availability of this 
amino acid is a determinant factor for growth capabil- 
ities of strain PB12.SA22 after 9 h of fermentation, al- 
though the strain could use additional (non-specified) 
carbohydrates supplied by YE (17.5%) [36] or other 
amino acids as possible carbon sources. 

Based on our GTA data, differential upregulation of 
genes encoding glycine transport components (operon 
proVWX) and the genes involved in the synthesis of 
serine from 3-P-glycerate and its subsequent transform- 
ation to glycine by upregulation oiglyA and serA, no differ- 
ential expression of serB and serC genes (Additional file 1) 
and the upregulation encoding genes for the GCV sug- 
gest an important role for these two amino acids 
in strain PB12.SA22 because they are most likely used 
additionally as protein components to yield C 2 units, 
which provide precursors for diverse important cellular 
process; indeed, it was reported that during growth on 
glucose, E. coli employs ~ 15% of the carbon assimilated 
in serine and its metabolites, such as the synthesis of 
cysteine, phospholipids, and glycine, whereas glycine is 
used in the synthesis of purines and heme-containing 
compounds. Ci units derived from serine and glycine 
are used in the synthesis of purines, histidine, thymine, 
pantothenate, and methionine and in the formylation of 
the aminoacylated initiator fMet-tRNA used to start 
translation in E. coli. In turn, these compounds are in- 
volved in other essential cellular process demanding the 
supply of Ci units [41,42], suggesting that particularly 
in the EXP/STA1 and EXP/STA2 comparisons Ci units 
derived from serine and glycine could be channeled for 
the biosynthesis of diverse cellular components used 
mainly for cellular maintenance of strain PB12.SA22 
(Figure 3). 

Interestingly, the genes encoding the complete biosyn- 
thetic pathway for arginine from glutamate ArgA, ArgB, 
ArgC, ArgD, ArgE, ArtI, ArgG and ArgH {argA, argB, 
argC, argD, argE, argF, argG and argH, respectively) and 
from bicarbonate to carbamoyl-phosphate (carAB) [17] 
were found to be upregulated in the STA1/STA2 com- 
parison, whereas they were downregulated during the 
EXP phase (Additional file 1), suggesting an important 
role of this amino acid during the STA phase distinct 
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Table 5 Average expression data of upregulated genes during STA1/STA2 comparison based on resultant FDR value 
from the RP analysis equal or less to 0.05 in compared growth phases during SA production cultures 



Gene 


Locus 


Protein name 3 


Cellular function GO b 


Average 
expression 


argA 


b2818 


Acetylglutamate synthase 


Amino acid metabolism 


87.9837 


argB 


b3959 


N-acetylglutamate kinase 


Amino acid metabolism 


45.8143 


orgC 


b3958 


N-acetylglutamylphosphate reductase 


Amino acid metabolism 


53.4932 


argD 


b3359 


N-succinyldiaminopimelate-aminotransferase/acetylornithine transaminase 


Amino acid metabolism 


21.9053 


argE 


b3957 


Acetylornithine deacetylase 


Amino acid metabolism 


18.7485 


orgG 


b3172 


Argininosuccinate synthase 


Amino acid metabolism 


65.2192 


argti 


b3960 


Argininosuccinate lyase 


Amino acid metabolism 


23.7859 


argl 


b4254 


Ornithine carbamoyltransferase chain I 


Amino acid metabolism 


33.2672 


corA 


b0032 


Polypeptide: CarA 


Amino acid metabolism 


107.7785 


corB 


b0033 


Polypeptide: CarB 


Amino acid metabolism 


23.9781 


ilvGJ 


b3767 


Acetolactate synthase II, large subunit, N-ter fragment 


Amino acid metabolism 


75.4038 


ilvG_2 


b3768 


Acetolactate synthase II, large subunit, N-ter fragment 


Amino acid metabolism 


24.6648 


IlvM 


b3769 


Polypeptide: IlvM 


Amino acid metabolism 


16.9332 


metE 


b3829 


Cobalamin-independent homocysteine transmethylase 


Amino acid metabolism 


59.0097 


serA 


b2913 


a-Ketoglutarate reductase/D-3-phosphoglycerate dehydrogenase 


Amino acid metabolism 


15.4266 


yboS 


b0485 


Glutaminase 


Amino acid metabolism 


27.4637 


fimC 


b4316 


Periplasmic chaperone, required for type 1 fimbriae 


Cell envelope 


20.5382 


ydiO 


b1697 


Putative subunit of YdiQ-YdiR flavoprotein 


Electron carrier activity 


32.4583 


ydiR 


b1698 


Putative subunit of YdiQ-YdiR flavoprotein 


Electron carrier activity 


32.3715 


ydiT 


b1700 


Putative ferredoxin 


Electron carrier activity 


28.0548 


ornB 


b2253 


UDP-L-Ara40 C-4' transaminase 


Lipid metabolism 


20.7697 


arnC 


b2254 


Undecaprenyl phosphate-L-Ara4FN transferase 


Lipid metabolism 


33.1058 


yhiD 


b3508 


Predicted Mg(2+) transport ATPase 


Lipid metabolism 


52.2858 


norY 


b1467 


Nitrate reductase Z, (3 subunit 


Other metabolic process 


18.7104 


narZ 


b1468 


Nitrate reductase Z, a subunit 


Other metabolic process 


13.4094 


yciE 


b1257 


Conserved protein 


Other metabolic process 


18.9020 


ydiS 


b1699 


Putative flavoprotein 


Other metabolic process 


20.0646 


ygdl 


b2809 


Putative lipoprotein 


Out of OGs 


15.8477 


sip 


b3506 


Starvation lipoprotein 


Out of OGs 


16.3895 


yboY 


b0453 


Predicted OM lipoprotein 


Out of OGs 


31.4852 


yhjR 


b3535 


Conserved protein 


Out of OGs 


17.5298 


yhiE 


b3512 


GadE DNA-binding transcriptional activator 


Regulatory 


224.6027 


yhiF 


b3507 


Polypeptide: predicted DNA-binding transcriptional regulator 


Regulatory 


33.7016 


yhiW 


b3515 


GadW DNA-binding transcriptional dual regulator 


Regulatory 


33.4367 


ykfE 


b0220 


Inhibitor of vertebrate C-type lysozyme 


Regulatory 


40.8144 


godB 


b1493 


Glutamate decarboxylase B 


Stress response 


94.5462 


hdeA 


b3510 


HdeA dimer, inactive form of acid-resistance protein 


Stress response 


125.3695 


hdeB 


b3509 


Acid stress chaperone 


Stress response 


104.9129 


hdeD 


b3511 


Acid-resistance membrane protein 


Stress response 


140.6337 


katE 


b1732 


Heme d synthase/hydroperoxidase 


Stress response 


15.6533 


cysD 


b2752 


Polypeptide: CysD 


Sulfur metabolism 


34.5193 


cysH 


b2762 


3'-Phospho-adenylylsulfate reductase 


Sulfur metabolism 


15.6801 


hisP 


b2306 


Lysine/arginine/ornithine ABC transporter/histidine ABC transporter - ATP binding subunit 


Transport 


1 7.8485 
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Table 5 Average expression data of upregulated genes during STA1/STA2 comparison based on resultant FDR value 
from the RP analysis equal or less to 0.05 in compared growth phases during SA production cultures (Continued) 



potG 


b0855 


Putrescine ABC transporter - ATP binding subunit 


Transport 


18.0562 


artJ 


b0860 


Arginine ABC transporter - periplasmic binding protein 


Transport 


120.3908 


gadC 


b1492 


Glutamic acid: 4-aminobutyrate antiporter 


Transport 


28.9528 


hisM 


b2307 


Lysine/arginine/ornithine ABC transporter/histidine ABC transporter - membrane subunit 


Transport 


44.5442 


norU 


b1469 


NarU MFS nitrate/nitrite antiporter 


Transport 


62.1767 


ompC 


b2215 


OMP C 


Transport 


37.5373 


yggB 


b2924 


Mechano sensitive channel MscS 


Transport 


22.3863 



Retrieved from EcoCyc database, b Biological functions were assigned according to the EcoCyc database for E. coli strain MG1655 and broadly grouped on the 
basis of GO terms [17]. 



from its role as a protein component. Arginine is also 
available in the supernatant broth because it is also sup- 
plied by YE (3.03%) [36]. This amino acid is transported 
by the arginine ABC transporter ArtPMQJI (argP, argM, 
argQ, argj and argl, respectively); however, among these 
genes, artl encodes the periplasmic binding component 
and is the unique found upregulated gene of this trans- 
porter, whereas the remaining genes were found without 
changes in their expression, suggesting that in addition 
to the upregulation of genes encoding the complete bio- 
synthetic pathway for arginine from glutamate and bicar- 
bonate, this amino acid is also transported from the 
fermentation broth during the STA phase. Arginine can 
be degraded by two pathways, the arginine succinyl- 
transferase (AST) pathway yielding succinate, which en- 
ters to TCA, and the arginine decarboxylase/agmatinase 
pathway, which also yields succinate. This last pathway 
yields the intermediate putrescin, which connects argin- 
ine degradation with the biosynthetic pathways of the 
polyamines putrescine, spermidine, cadaverine and ami- 
nopropylcadaverine [17]. Interestingly, unlike E. coli K-12, 



wild-type E. coli strains are unable to use L-arginine as 
a carbon source, but they can use this amino acid as a 
nitrogen source [47,48]. Based on this evidence, our tran- 
scriptomic data suggest that arginine could be used by 
strain PB12.SA22 as a source of succinate to recycle C 
through TCA or as a source for polyamines during the 
STA phase for SA production. However, data mining of 
the RP method analysis for microarray data comparison 
between STA1/STA2 showed that none of the genes of 
the arginine AST degradation (operon astCADBE) and ar- 
ginine decarboxylase/agmatinase (adiA, speB, speF, puuAP, 
puuDRCBE and sad) pathways were significantly upregu- 
lated (Additional file 1). Polyamines are reported as neces- 
sary for cell growth because they are major polycations in 
cells, together with Ca 2+ and Mg 2+ . Polyamines and Mg 2+ , 
which are present in higher free concentrations than Ca 2+ , 
can bind to intracellular polyanions such as nucleic acids 
and ATP to modulate their function [49] . Polyamines have 
also been associated with stabilizing membranes and 
stimulating several enzymes [17], and furthermore, spermi- 
dine was found to donate a portion of its molecule for the 




Figure 3 Proposed metabolic and physiologic traits of E. coli PB12.SA22.based on relevant upregulated genes during EXP and STA 

GTA. OM, IM, Glycolysis, PPP, TCA and SA pathway are illustrated schematically. Arrows indicate the direction of transport or reaction; dotted 

arrows indicate suggested or unknown mechanisms. 
\ ) 
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enzymatic biosynthesis of hypusine, a unique amino acid 
that plays a crucial role in cell proliferation [50]. Because 
strain PB12.SA22 ceased EXP growth by hour 9 of cultiva- 
tion and entered a pseudo STA phase until the end of the 
fermentation, it is unlikely that arginine could be used as a 
polyamine precursor, suggesting that it could be used to 
supply succinate to TCA, through the gluconeogenic cap- 
abilities developed by strain PB12 [18,51], to make PEP 
available for use in SA production, particularly during the 
late STA phase in the absence of glucose (Figure 3). 

Possible associations between amino acid transport, 
catabolism and pH stress response 

Amino acids such as aspartate, serine and tryptophan 
could be used potentially by E. coli as carbon sources 
under certain growth conditions. Catabolism of aspartate 
by aspartate ammonia-lyase, AspA (aspA) yields Fum + 
ammonia + H + ; catabolism of L-serine by serine deami- 
nase, SdaB (sdaB) yields PYR + ammonia + H + ; whereas 
the catabolism of tryptophan catalyzed by tryptophanase 
(tnaA) yields PYR + ammonia + indole + H + . In turn, fu- 
marate and PYR could fuel TCA or increase intracellular 
PYR availability, respectively [17]. However, amino acid 
deamination is a metabolic process displayed by E. coli to 
contend with alkaline stress that supplies carbon for bac- 
terial growth under this condition, particularly when 
growing in complex peptide-rich media. As cell density in- 
creases, extracellular pH alkalinizes and released ammonia 
is deprotonated and volatilized, whereas the C skeleton is 
channeled into acids [52-54]. Based on our transcriptomic 
data showing the upregulation of genes such as tnaA, 
aspA and sdaB during the EXP and STA phases, we 
propose that intracellular availability of amino acids such 
as aspartate, cysteine and tryptophan, either by importing 
them from the extracellular environment (as they are pro- 
vided by YE) or by interconversion from other sources, 
e.g., cysteine from O-acetyl-L-serine by CysK (cysK) [17] 
(also upregulated during the EXP/STA2 comparison), 
have an important role in a possible cellular response of E. 
coli PB12.SA22 during SA production to extracellular 
alkalinization because upregulation of tnaA, asp A and 
sdaB suggests that their respective enzymes deaminate 
aspartate, cysteine and tryptophan to PYR or Fum. In 
addition to the upregulation of these genes, a previous re- 
port indicates that as a consequence of cellular exposure 
to neutral pH conditions, E. coli 3110 showed a higher 
expression of 93 genes (grouped as the neutral high 
expressed cluster showing highest expression at pH 7.0 
and lower expression at both pH extremes), including 
fecA, fecB, fecE, lamB, UvK, malE, malF and malK [55]. 
Interestingly, all these genes were found to be upregulated 
in strain PB12.SA22 in the EXP/STA1 comparison, with 
exception of malE, which was also upregulated in the 
STA1/STA2 comparison. 



Several genes encoding proteins involved in the AR in 
E. coli were found to be upregulated during the EXP and 
STA phases, including the chaperone-based AR proteins 
HdeA and HdeB, which belong to the hdeAB acid stress 
operon involved in the periplasmic acid stress response 
that prevent periplasmic-protein aggregation at acidic 
pH; the IM protein HdeD (hdeD), which is required for 
acid resistance [56], was found to be upregulated during 
the STA1/STA2 comparison, and GadB and GadC 
(gadBC operon) were upregulated in the STA1/STA2 
phases and are components of the Gad system (the glu- 
tamate decarboxylase system), the major AR system in 
E. coli under extreme acidic conditions [17]. GadB (to- 
gether with GadA) are pyridoxal 5 '-phosphate ^In- 
dependent enzymes that convert glutamate to y-amino 
butyric acid (GABA) and carbon dioxide in a reaction 
that consumes a cytoplasmic proton. GABA is trans- 
ported out of the cell by the IM antiporter GadC in ex- 
change for more extracellular glutamate [57], whereas 
GadW is a transcriptional dual regulator (yhiW or 
gadW) [17,57,58]. Several regulatory genes of the Gad 
system (gadA, gadE, gadW and gadX), as well as the 
chaperone-related genes hdeABD, were placed in a 
unique cluster termed the acid fitness island [59], but 
interestingly, more genes in this island were found to 
participate in acid resistance, including the OM lipopro- 
tein Sip (sip) and the transcriptional regulator YhiF 
(yhiF) (both genes were upregulated in the STA1/STA2 
comparison in our transcriptomic study), which are re- 
quired to protect cells against excreted toxic metabolites, 
including the accumulated anions of dissociated weak 
acids after growth at low pH such as lactate, succinate, 
and formate [57]. Additionally, GadE serves as the global 
transcriptional activator for many genes involved in 
stress response, glutamante biosynthesis, and in the bio- 
synthesis of membrane components. Auto-induction of 
GadE requires the alternative sigma s factor responsible 
for the transcription of many genes in the STA phase 
and another unidentified factor, whereas P2 and P3 of 
gadE are activated by GadX (or YhiX) and GadW (or 
YhiW) during the STA phase [57]. 

The simultaneous upregulation of genes encoding pro- 
teins involved in acid pH stress response with the upregu- 
lation of catabolic amino acid pathway to revert an 
apparent alkalinization has been reported to have a com- 
mon connection, the protein GadB [52]. As was discussed 
above, high extracellular pH induces transcription of 
several genes encoding enzymes generating ammonia, 
such as TnaA, CysK, SdaB and AsnA, by deaminating 
tryptophan, serine, cysteine and aspartate, respectively. Ex- 
pression of gadB is controlled by sigma s during the STA 
phase in minimal growth medium; however, induction of 
this gene has been observed during growth in complex 
medium between pH 5.5 and 8.0 [52,55,60]. Induction of 
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the RpoS dependent- oxidative AR system has been re- 
ported due to the presence of glutamate in YE. Glutamate 
and glutamine appear to activate a preformed RpoS- 
dependent system that is produced by entry into the STA 
phase [60]. According to this observation, simultaneous 
upregulation of the AR Gad system associated with the in- 
duction of genes involved in alkaline pH stress observed 
during SA production could be induced by the presence 
of glutamate in YE, which is transported by GadC. Add- 
itionally, upregulation of hdeAB during the EXP/STA2 
comparison can be explained by the inductor effect of glu- 
tamate because these genes are also activated by GadE 
and GadW [17]. 

It has been established that the OM envelope and the 
periplasmic space are cellular compartments exposed es- 
sentially to extracellular pH [52]. In bacterial strains 
used for metabolite production such as PB12.SA22, dur- 
ing SA production, it is not surprising that the periplas- 
mic side of the IM, the periplasmic space and the inner 
side of the OM could be exposed to low acid stress be- 
cause SA is exported from the cytoplasm and then dif- 
fuses into the extracellular environment; conversely, the 
external side of the OM could be exposed to low alka- 
line stress as a consequence of the addition of NH 4 OH 
to maintain fermentor pH = 7.0. This supposition makes 
sense given the observed upregulation of genes involved 
in acid pH stress response, simultaneous with the upreg- 
ulation of the deaminating amino acid pathway to re- 
verse an apparent alkalinization during SA production 
by strain PB12.SA22 (Figure 3). 

Upregulation of genes involved in iron and sulfur 
transport and its metabolism 

Iron and sulfur seems to play an important physiological 
role during SA acid production by strain PB12.SA22 be- 
cause some genes encoding transporters for these ions 
were found to be upregulated. Among these, fecBE was 
found to be upregulated in the EXP/STA1 comparison, 
and fecA was observed in the EXP/STA2 comparison. 
The fee ABC operon genes are located in the fecABCDE 
operon and encode the ferric citrate OMP FecA (fecA), 
the IM ferric dicitrate ABC transporter (fecB), and the 
ATP binding component (fecE) [61]. Additionally, it was 
found to be upregulated in the EXP/STA1 comparison 
the OM receptor Cir (cirA), a TonB dependent iron- 
siderophore complex involved in iron uptake and regu- 
lated by both cellular iron content and growth [62]. Re- 
garding sulfur, the upregulation of the periplasmic 
binding protein (Sbp) sulfate/thiosulfate ANC trans- 
porter (sbp) in the EXP/STA1 comparison was observed, 
which is associated with the transport both of sulfate 
and thiosulfate that are used as sulfur sources [63]. Fi- 
nally, in the STA1/STA2 comparison, we also found that 
the genes narll, narZ and narY (part of the narllZYWV 



operon) were upregulated, which encode the nitrite/ni- 
trate transporter (narll) and the a and ft subunits of the 
IM located nitrite /nitrate reductase (narZ and narY, 
respectively). 

Iron is essential for the growth of E. coli, as is empha- 
sized by the variety of processes in which iron- 
containing proteins take part, including their structural 
association to proteins involved in electron transport or 
playing important roles, particularly for iron-sulfur pro- 
teins involved in amino acid and pyrimidine biosynthesis 
(glutamate synthase, dihydroo rotate dehydrogenase), and 
the TCA (aconitase, succinate dehydrogenase), as well as 
in electron transport (ferredoxin) and non-heme, non- 
iron-sulfur proteins required for DNA synthesis (ribonu- 
cleotide reductase), protection from superoxide radicals 
(superoxide dismutase), and interestingly, aromatic 
amino acid biosynthesis (DAHP synthase) [64]. Exogen- 
ous ferric citrate (supplied in the fermentation broth for 
SA production as ammonium iron (III) citrate), is trans- 
ported across the OM by FecA, and a signal is transmit- 
ted across the OM to the IM protein FecR, which 
transmits the signal across the IM, thereby activating 
(through the sigma 70 factor) the cytoplasmic family pro- 
tein Feci, which directs RNA polymerase to express the 
fecABCDE operon [17,61,65]. 

Upregulation of genes involved in nucleotide/nucleoside 
transport and biosynthesis 

During the EXP/STA1 comparison, upregulation was ob- 
served for the gene encoding OMP Tsx (tsx), a protein 
that has been proposed to function in vivo as a pore that 
specifically facilitates the permeation of nucleosides and 
deoxynucleosides across the OM due to its specificity for 
free nucleobases or monophosphate nucleosides [22,66]. 
Regarding nucleotide/nucleoside metabolism, were found 
upregulated genes involved in pyrimidine ribonucleoside 
salvage, rihA and upp encoding ribonucleoside hydrolase 
1 (pyrimidine-specific) and uracil phosphoribosyltransfer- 
ase, respectively, which catalyze the sequential conversion 
of uridine — > uracil — > UMP [17]. Regarding pyrimidine ri- 
bonucleotide de novo biosynthesis, was found upregulated 
the gene pyrB encoding the catalytic subunit of the aspar- 
tate carbamoyltransferase, an enzyme catalyzing the con- 
version of L-asparte + carbamoyl-P to N-carbamoyl-L- 
aspartate in reactions involved in the biosynthesis of 
UMP in this pathway [17]. Finally, regarding the aden- 
ine and adenosine salvage pathway, upregulation was 
detected for the gpt gene, encoding xanthine-guanine 
phosphoribosyltransferase, which catalyzes the conver- 
sion of hypoxanthine to IMP, a common precursor for 
the guanosine and adenosine nucleotides in the de novo 
biosynthetic pathway [17]. 

Xanthine-guanine phosphoribosyltransferase is also in- 
volved in the xanthine and xanthosine salvage pathway, 
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catalyzing the conversion of xanthine to XMP, which is 
then channeled to the purine nucleotide de novo biosyn- 
thesis pathway, and also catalyzes the transformation of 
guanine to GMP in the guanine and guanosine salvage 
pathway. GMP is, in turn, channeled to the guanosine nu- 
cleotide de novo biosynthesis pathway [17]. Additionally, 
upregulation was observed for guaC, encoding the GMP 
reductase, which catalyzes the conversion of GMP to IMP 
by a reductive deamination, and the nrdA gene, encoding 
the a-subunit dimer of the ribonucleoside diphosphate re- 
ductase 1 enzyme, which catalyzes the conversion of nu- 
cleotides to deoxynucleotides, an essential step during 
DNA synthesis, including its role in the chromosome rep- 
lication and repair processes [17]. Data mining of the 
microarray expression comparison between the EXP/ 
STA1 phases showed no significant expression differences 
based on the FDR value < 0.05 for other genes involved in 
salvage pathway of pyrimidine ribonucleotides (pyrH, ndk, 
pyrG, cmk) or for the pyrimidine de novo biosynthesis 
pathway (pyrC, pyrD, pyrE, pyrF, ndk or pyrG) (Additional 
file 1). Interestingly, nrdA was found to be upregulated in 
the late STA phase, but no other genes involved in the 
purine and pyrimidine pathways discussed above were 
found to be significantly upregulated in the STA1/STA2 
comparison. 

All genes encoding enzymes related to the so-called 
nucleotide salvage pathways are involved in transforma- 
tions of purine nucleotides derived from exogenous 
sources, whereas in the de novo pathways, they are syn- 
thesized from simpler precursors [67]. Our transcrip- 
tomic data, showed the upregulation of few genes 
encoding enzymes catalyzing reactions both in the de 
novo and salvage pathways involved in the synthesis 
of the purines IMP, XMP, GMP and the pyrimidine nu- 
cleotide UMP. Upregulation of the above described 
genes ensures the availability of rNTPs and dNTPs for 
RNA and DNA synthesis, respectively, from both purine 
(IMP, XMP, GMP) and pyrimidine (UMP) essential nu- 
cleotides during the growth stage of strain PB12.SA22 
(Figure 3). 

Modification of the cellular membrane during the STA 
phase 

In the STA1/STA2 comparison, important differences 
were observed in the upregulation of several genes pos- 
sibly involved in inner and OM structure modification 
and the response to diverse environmental processes 
with respect to the EXP/STA1 comparison. These in- 
clude the arnBC genes encoding an undecaprenyl trans- 
ferase involved in the modification of lipid A phosphates 
with 4-amino-4-deoxy-L-arabinose (L-Ara4N), which 
modifies the lipid composition of the outer face of the 
periplasmic membrane [68]; yciE, encoding YciE, a con- 
served protein of unknown function, which has been 



observed to be induced under osmotic stress imposed by 
NaCl in both aerobic and anaerobic conditions [69]; 
ygdl, encoding an IM putative lipoprotein [17]; sip, 
encoding a STA phase lipoprotein (Sip) that has been 
proposed to take part in acid resistance because its ex- 
pression was observed when cells were grown in pH 5.5 
to 4.5 under conditions known to induce glutamate 
dependent acid resistance compared to pH 7.4, and also 
associated with YhiF to protect the cell against toxic me- 
tabolites [54]. Expression of sip was previously observed 
to increase 3-5 -fold in the STA carbon starvation and 
was found to form homo-oligomer complexes tethered 
to the OM [70]. Other genes found to be upregulated 
were ydhl (or yhhE), a putative ATPase transporter in- 
volved in high cell dependent acid resistance [71]; and 
the ompC gene, encoding OmpC. Although the possible 
role of OmpC in glucose transport during EXP growth 
was discussed above, expression of ompC and OmpC 
protein level have been demonstrated to also be influ- 
enced by a wide variety of environmental conditions in- 
cluding pH, osmolarity, temperature, concentration of 
certain toxins, and growth phase [72]. OmpC and OmpF 
are reported as the major constituents of the OM in 
E. colU accounting for approximately 2% of the total pro- 
tein content of the cell [73]. The role of the OMPs in 
the modification of the OM has been previously re- 
ported, indicating that the OmpC/OmpT:OmpA ratio 
increased in total membrane protein analysis of an 
evolved E. coli resistance to increasing isobutanol con- 
centrations correlating with upregulation of these genes, 
resulting in changes in the OM structure, accomplished 
with modification in membrane composition and pep- 
tidoglycan structure [74]. Interestingly, in support of our 
proposition regarding the possible role of those genes 
discussed in this section, particularly during the STA 
phase, cell resistance to lysis during the total RNA ex- 
traction procedure showed that PB12.SA22 modified its 
cellular surface properties because the protocol origin- 
ally developed for successful total RNA extraction in the 
mid EXP phase (5 h) was modified for the extraction in 
the STA phases, thereby enhancing the cellular lysis step 
to yield a high amount of total RNA suitable for micro- 
array analysis. 

Sigma factor regulatory networks controlling the 
expression profiles in the EXP and STA phases 

We constructed each regulatory network of the sigma 
factors controlling upregulated genes observed in EXP/ 
STA1, EXP/STA2 and STA1/STA2 comparisons. The 
sigma regulatory interactions were retrieved from the 
comparisons between the information stored in the Eco- 
cyc and Regulon DB databases [17,45]. The comparisons 
show that the sigma 70 (RpoD) was the main sigma factor 
targeting promoters for genes encoding amino acid 
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metabolism, ATP biosynthetic/catabolic process, trans- 
port, electron carrier activity, nucleotide/nucleoside me- 
tabolism, cell envelope and stress response (Figure 4). 
However, although RpoD is the primary sigma factor 
during EXP growth, targeting a wide range of promoters 
that are essential for normal growth [75] as is shown in 
Figure 4, transcription of several upregulated genes were 
possibly controlled by sigma 70 and at least one or two 
more sigma factors. However, the transcription of other 
sets of genes are possibly controlled only by alternative 
sigma factors, and interestingly, 17%, 16% and 18% of 
significantly upregulated genes during the EXP/STA1, 
EXP/STA2 and STA1/STA2 comparisons, respectively, 
are controlled by unspecified sigma factors. 

Alternative sigma factors possibly controlling tran- 
scription of several genes simultaneously to sigma 70 dur- 
ing EXP/STA1 comparison are sigma s , sigma 32 and 
sigma 19 . Interestingly, the fecA, fecB and fecE genes (en- 
coding the ferric citrate OMP, the periplasmic binding 
protein and the ATP binding subunit of the dicitrate 
ABC transporter, respectively), are located in an operon 
with a promoter targeted by sigma 19 controlling the ex- 
pression of the fecAB CD operon genes, whereas the ex- 
pression of the contiguous gene in the operon, fecE> is 
controlled by sigma 70 and sigma 19 , thus suggesting, as 
discussed above, an important role for iron during EXP 
growth of strain PB12.SA22. During EXP/STA2 com- 
parison, sigma s and sigma 70 were found to be control- 
ling the expression of six genes involved in transport 
(membrane subunit of the glycine betaine/proline ABC 
transporter encoded in the operon proVWX), acid stress 
response (hdeAB operon) and DNA binding transcrip- 
tional activators (yhiE, yhiW genes), respectively [17]. 
Data retrieved from the EcoCyc and Regulon DB data- 
bases suggested that during the STA phase (the STA1/ 
STA2 comparison), sigma 70 was most likely associated 
with the transcription of 35 upregulated genes; however, 
in this comparison, it was proposed that the sigma 



factors 54 , s , 32 and 24 possibly simultaneously participated 
with sigma 70 RpoD in the co-transcription of 12 genes, in- 
cluding those involved in acid stress response (gadBC, 
hdeAB and yhiD), nitrite/nitrate transport (yciE), osmotic 
stress response (yciE, katE) and DNA binding transcrip- 
tional activators (yhiE, yhiW genes) [17] (Figure 4). 

Sigma 70 accounts for 60-95% of the total pool of cellular 
sigma factors during normal EXP growth [17]. Changes to 
typical growth conditions, such as heat shock, acid stress 
or growth into the STA even in rich broth, lead to the re- 
placement of RpoD with other sigma factors such as RpoS, 
which is considered to be the master regulator of the gen- 
eral stress response in E. coli. RpoS is practically absent in 
rapidly growing cells but is strongly induced during entry 
into the STA phase and/or many other stress conditions 
and is essential for the expression of multiple stress resis- 
tances [76]. When wild-type E. coli is grown in glucose- 
limiting conditions and, interestingly, in PTS~ strains such 
as PB12.SA22 parental strain PB12, transcription of sev- 
eral central metabolism and especially glycolytic genes 
turn on under the control of RpoS, which has been pro- 
posed as a second vegetative sigma factor with major im- 
pact not only on stress tolerance but on the entire cell 
physiology under non-optimal growth conditions [77,78]. 
However, the comparison of our transcriptomic data 
showed that the significantly upregulated genes across 
EXP/STA1, EXP/STA2 and STA1/STA2 that are con- 
trolled by RpoD average 71%, and interestingly, in the 
comparison between EXP/STA2 and STA1/STA2, the 
sigma s factor was involved in the simultaneous upregula- 
tion of 22% of those genes controlled by sigma 70 . 

Data mining of differentially expressed genes during 
all performed comparisons showed that all genes encod- 
ing sigma factors were not significantly up- or downreg- 
ulated based on the FDR value < 0.05 (Additional file 1). 
This result suggests a relevant role for those positive or 
dual transcriptional factors targeting transcription initi- 
ation sites of upregulated genes. 




Figure 4 Proposed transcriptional regulatory network of the sigma factors in E. coli PB12.SA22 in EXP/STA1, EXP/STA2 and STA1/STA2 
comparisons. U, unspecified sigma factor. 



Cortes-Tolalpa et al. Microbiol Cell Factories 2014, 13:28 
http://www.microbialcellfactories.eom/content/13/1/28 



Page 16 of 20 



Concluding remarks 

During batch cultures for SA production using complex 
fermentation broth (including 25 g/L of glucose and 
15 g/L of YE), strain PB12.SA22 ceased EXP growth 
even in the presence of a high amount of residual glu- 
cose, indicating that growth was not associated with glu- 
cose consumption, depending to some extent of the 
availability of nutrients supplied by YE. This hypothesis 
was supported by the observation that EXP growth ces- 
sation is associated with the total consumption of avail- 
able tryptophan in the supernatant, entering a STA-like 
phase where a brief increment in biomass was observed. 
During this stage, cells consumed residual glucose in the 
supernatant, suggesting that biomass produced during 
the EXP phase entered a resting cell-like condition pro- 
ducing SA. Because no changes were observed in the 
regulation of genes involved in CCM and the SA 
pathway between the EXP and STA phases, it is pos- 
sible to propose that this strain transports and catab- 
olizes extracellular glucose mainly for SA production 
in the STA phase until 26 h of cultivation, where 
glucose was completely consumed. Those metabolic 
capabilities regarding glucose catabolism and synthesis 
of precursors PEP and E4P for SA production could 
be a consequence of the permanent scavenging condi- 
tion proposed previously for parental strain PB12, as 
a consequence of the inactivation of the PTS operon 
and diverse genetic changes developed during an evo- 
lutive adaptation process from which the PB12 strain 
was obtained [18,79]. 

GTA of E. coli strain PB12.SA22 during SA production 
between the mid EXP phase (5 h), the early STA phase 
(9 h) and the late STA phase (44 h) showed no signifi- 
cant differences in absolute expression in genes encod- 
ing enzymes of CCM, glycolysis, the PPP, TCA, the 
glyoxilate shunt, acetate metabolism and gluconeogenic 
enzymes, such as those in the SA pathway. However, im- 
portant differences were observed in the upregulation of 
genes encoding proteins involved in sugar, amino acid, 
nucleotide/nucleoside and iron and sulfur transport; 
metabolic processes, particularly amino acid catabolism 
and biosynthesis; nucleotide/nucleoside salvage; acid 
stress response; and the modification of the cell mem- 
brane in the EXP/STA1, EXP/STA2 and STA1/STA2 
comparisons. Figure 3, shows a proposed metabolic 
model based on relevant upregulated genes observed 
during comparisons in strain PB12.SA22 during SA pro- 
duction in batch-fermentor cultures grown in a complex 
broth. During the EXP/STA1 comparison, strain PB12. 
SA22 has an important transport activity through OMPs 
and IM proteins that possibly transport glucose, sulfur, 
iron, amino acids, peptides and nucleotide/nucleosides 
whereas amino acid interconversion and catabolism 
process were maintained. In the EXP/STA2 comparison, 



an important reduction in transport activity was ob- 
served particularly for some amino acids and possibly 
peptides supplied by the YE, while important amino acid 
interconversion and catabolism processes, such as nu- 
cleotide/nucleoside biosynthesis, were maintained. In 
both comparisons, it is important to highlight the catab- 
olism of amino acids such as serine, tryptophan and as- 
paragine (EXP/STA1) and serine and tryptophan (EXP/ 
STA2), which results in PYR or Fum, and particularly 
the possible catabolism of glycine through the GCV, sug- 
gesting the possibility that PYR and Fum could be chan- 
neled to the CCM, whereas the catabolism of cysteine 
could result in the availability of Q units for the biosyn- 
thesis of cellular blocks. Finally, during the STA1/STA2 
comparison, only the transport of arginine, glutamic 
acid, and nitrite/nitrate and the biosynthesis of arginine 
from glutamic acid were observed. The availability of ar- 
ginine and its possible catabolism could yield succinate, 
which could be channeled to the TCA. As was discussed, 
glucose is completely consumed at 26 h of cultivation, 
but SA accumulation continued until the end of fermen- 
tation. The results suggest how strain PB12.SA22 con- 
tinues producing SA during the late STA phase. Even if 
the catabolism of arginine could fuel carbon into TCA 
and, possibly, the gluconeogenic capabilities of this 
strain could supply PYR and PEP, the availability of E4P 
becomes a limiting resource for DAPH synthesis and the 
subsequent accumulation of SA. Another important 
group of genes upregulated in the STA1/STA2 com- 
parison are possibly involved in the cellular response 
to pH stress genes (gadB, gadC, hdeA, hdeB) and 
several genes encoding proteins possibly involved in 
both outer and IM modifications as a response to en- 
vironmental conditions imposed on the cell at the 
end of the fermentation. 

Finally, the network reconstruction based on different 
sigma factors controlling the expression of upregulated 
genes of each condition showed that even sigma s has 
been proposed to be the master regulator of the STA, 
and sigma 70 also plays an important role in controlling 
some upregulated genes with no alternative reported 
sigma factor or any sharing its regulatory region with a 
sigma s promoter or an alternative sigma factor. This 
sharing is an important property of regulatory networks, 
which gives them plasticity to adapt to different environ- 
mental conditions. 

Methods 

Bacterial strain and growth conditions 

The SA-producing strain PB12.SA22 (PTS" AaroK AaroL 
pTOPOaroBaroE, pILBaroG^tktA) [10] was shake-flask 
cultured in 125 mL baffled flasks containing 10 mL of Luria- 
Bertani broth supplemented with 30^g/mL of kanamycin, 
15^g/mL of gentamycin, 20^g/mL of chloramphenicol, 
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and 30^g/mL of tetracycline. Cultures were incubated 
overnight in a shaker (New Brunswick Scientific, Edison, 
USA) at 37°C, 300 rpm. An aliquot of 150 yL of this cul- 
ture was used to inoculate a 250 mL baffled flask with 
50 mL of fermentation medium containing 25 g/L of glu- 
cose and 15 g/L of YE [8,10]. Biomass concentration was 
determined, and calculations were performed to adjust 
inoculum size at OD 600 nm = 0.35 and to inoculate batch 
cultures (by triplicate) in an Applikon autoclavable glass 
Bio Reactor one L fermentor with a 500 mL working vol- 
ume of fermentation medium supplemented with required 
antibiotics. Fermentor vessels were connected to an 
Applikon ADI 1010 BioController and ADI 1025 control- 
lers to monitor temperature, pH, impeller speed and DO. 
Batch fermentations were run for 50 h at 37°C, pH 7.0 
(maintained by addition of 3.0% of NH 4 OH). Impeller 
speed of no less than 500 rpm was used to maintain DO 
levels at 20% air saturation [80]. Gene-carrying plasmids 
were induced by adding 0.1 mM of isopropyl /2-D-l- 
thiogalactopyranoside. 

Analytical procedures 

Biomass concentrations were monitored each hour dur- 
ing the first eight h of cultivation. After this point, con- 
centrations were monitored every six h until the end of 
the fermentation (50 h). Triplicated samples (1.5 mL) 
were withdrawn from each reactor and OD 600 nm was 
determined spectrophotometrically (Beckman DUR-70 
Spectrophotometer). Supernatant samples for the deter- 
mination of SA, DHP, DHS, and GA were prepared by 
centrifuging one mL of fermented broth at 12,000 rpm 
for 1 min. To remove any residual cell in the super- 
natant, the centrifuged samples were filtered through 
0.45 nylon membranes. SA, DHS, and GA concentra- 
tions were determined by HPLC as described previously 
[10], whereas DAHP concentrations were determined by 
the thiobarbituric acid assay [81]. Because this method 
does not distinguish between DAHP and DAH, in this 
work, DAHP levels corresponded to the sum of both 
compounds [82]. Glucose concentration was assessed 
by a biochemical analyzer (YSI 2700 Select). Aromatic 
amino acids phenylalanine, tyrosine and tryptophan, 
present in culture supernatants, were quantified using a 
Phenomenex Synergy Hydro RP18 column (150 by 
4.6 mm; 4/mi) attached to an Agilent 1100 HPLC sys- 
tem (Agilent Technologies). Running conditions were as 
follows: mobile phase, 0.2% of trifluoroacetic acid 40% 
of methanol; flow, 0.5 mL/min. Detection was per- 
formed by photodiode array at 280 nm [82]. 

The ft and the qs were calculated during the EXP 
growth phase, qs was calculated as the differential 
change in S, with time (t) normalized to the biomass 
concentration (#s = ^). A predetermined correlation 



factor (1 OD 600 nm corresponded to 0.37 g/L of DW) 
was used to transform OD 600 nm values into cell concen- 
trations for qs and \a calculation. 

TACY determinations were based on the combined 
molar yields of DAHP, DHS, SA, and GA [10]. Results 
were reported as the average of triplicated samples with 
their associated standard deviation. 

RNA extraction procedures 

One aliquot of 5 mL from each of three fermentor cul- 
tures performed were collected at the middle of the EXP 
growth phase (5 h, D0 60 o nm = 7.00), the early STA state 
(9 h, D0 60 o nm = 13.00) and the late STA state (44 h, 
D0 60 o nm =15.2) (resulting in triplicate samples for each 
time), and transferred to a 15-mL Falcon tube contain- 
ing 500 uL of RNA later solution (Ambion). Total RNA 
extraction was performed using a modification of the 
previously reported hot phenol-based procedure. Previ- 
ous methodology was used successfully for total RNA 
extraction from cultures of E. coli JM101, and its PTS~- 
derivative PB11 and PB12 strains, grown in batch- 
fermentor cultures in minimal M9 medium and collected 
at D0 60 o nm = 1.0 [18]; however, direct application of 
this methodology to biomass collected from cultures of 
strain PB12.SA22 for SA production in a complex 
medium resulted in poor total RNA yield, particularly in 
samples collected during STA. To obtain total RNA of 
high quality (non-degraded and with high absorbance 
DO 2 60/280 and DO230/260 ratios) and concentration re- 
quired to perform successful GTA, the original RNA ex- 
traction procedure was modified and optimized. The 
integrity of extracted total RNA was evaluated by agar- 
ose gel electrophoresis, and quality was evaluated by 
determination of concentration and the DO 2 60/280 and 
DO260/230 ratios in a Nanodrop-2000c spectrophotom- 
eter (Thermo Scientific). Samples were adjusted to a 
final concentration of total RNA = 1 fig/pL. Aliquots 
of 20 \aL of total RNA were shipped to the Precision 
Biomarker company for microarray experiments (www. 
precisionbiomarker.com). 

Microarray design, experiments and data analysis 

Total RNA extracted from each samples withdrawn at 5, 
9 and 44 h from triplicated fermentor cultures, were en- 
tirely processed by Precision Biomarker. Sample prepar- 
ation and processing included cDNA synthesis, cDNA 
fragmentation, and the preparation of the hybridization 
mixture, hybridization with the Affymetrix GeneChip® 
E. coli Genome 2.0 (Affymetrix), washing, staining, and 
scanning of the microarrays, which were all performed 
according to the manufacturers' procedures. Raw data 
acquisition and gene expression data analysis were per- 
formed by Precision Biomarker using background cor- 
rection, normalization and summarization. The resulting 
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databases corresponding to the analyzed samples were 
analyzed in our laboratory by retrieving expression data 
from the array of E. coli MG1655, including in the Affy- 
metrix GeneChip® E. coli Genome 2.0. Average data of 
triplicates from each growth phase sampled were used 
to calculate the relative expression value between EXP 
relative to STA1 (EXP/STA1 ( 2 a( EXP STA1 )), STA1 relative to 
STA2 (STA1/STA2 (2a( sta1STA2 )) an d EXP relative to STA2 
(2A (exp-sta2) ) The resulting databases were formatted using 
a Perl language program from http://www.gla.ac.uk/schools/ 
computing/?CFID=21255484&CFTOKEN=17693561 for 
analysis by the rank product method to identify differentially 
expressed (upregulated or downregulated) genes [15,16]. 

Additional file 



Additional file 1: Differentially upregulated and downregulated 
genes during EXP/STA1, EXP/STA2 and STA1/STA2 comparisons 
based on the application of the RP defined by a FDR value < 0.05. 
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